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In General Relativity without a cosmological constant a non-positive contribution from the space- 
time geometry to Raychaudhuri equation is found provided that particular energy conditions are 
assumed and regardless the considered solution of the Einstein's equations. This fact is usually 
interpreted as a manifestation of the attractive character of gravity. Nevertheless, a positive contri- 
bution to Raychaudhuri equation from space-time geometry should occur since this is the case in an 
accelerated expanding Robertson- Walker model for congruences followed by fundamental observers. 
Modified gravity theories provide the possibility of a positive contribution although the standard 
energy conditions are assumed. We address this important issue in the context of f(R) theories, 
deriving explicit upper bounds for the contribution of space-time geometry to the Raychaudhuri 
equation. Then, we examine the parameter constraints for some paradigmatic f(R) models in order 
to ensure a positive contribution to this equation. 



PACS numbers: 04.50.Kd, 95.36.+X, 98.80.-k 

I. INTRODUCTION 

Since observational evidence of the accelerated expan- 
sion of the Universe was discovered pQ, the cosmological 
evolution as predicted by General Relativity (GR) has 
been set in doubt. The reason is that a stress-energy 
tensor possessing strange features needs to be included 
in the Einstein's equations in order to account for this 
cosmic acceleration. This exotic cosmological fluid is usu- 
ally referred to as dark energy (DE). In its simple form 
it is given by a cosmological constant with equation of 
state pa = — Pa- However, instead of filling the Uni- 
verse with exotic fluids, a reasonable hypothesis consists 
in modifying the cosmological field equations assuming 
alternative geometrical theories to GR. This approach 
has received the name of modified gravity theories and 
has drawn some attention in the last years [2]. 

Some examples are Lovelock theories, whose field equa- 
tions are second-order differential equations in the met- 
ric [3J; Gauss-bonnet theories inspired in string theory 
that include a Gauss-Bonnet term in the Lagrangian 0] ; 
scalar-tensor theories [5] or vector-tensor theories [6], in 
which gravitational interaction is not only mediated by 
the standard spin-2 graviton but also by scalar or vector 
modes respectively; metric theories derived by extra di- 
mensional theories [7] ; supergravity models [8] , disformal 
theories jS], Lorentz violating and CPT breaking models 
of gravity TOJ ; or the so-called f(R) theories, in which 
our work will be focused. f(R) theories consist in mod- 
ifying the Einstein-Hilbert Lagrangian by adding an ar- 
bitrary function of the Ricci scalar R (for recent reviews 
see |11|). From this approach, the equations derived from 
the new action are to be expected as a refinement of the 
standard Einstein's equations able to reproduce the cor- 
rect predictions of GR while explaining the cosmic ac- 



celeration. These theories may have strong effects on 
small scales, but if some restrictions are imposed, they 
are able to reproduce the cosmological history while be- 
ing compatible with local gravity tests [12]. It is worth 
mentioning that the Einstein's equations with cosmologi- 
cal constant A are a particular case of these theories with 
f(R) = -2A. 

In fact, the problem with the accelerated expansion 
of the Universe follows immediately from the consid- 
eration of the Friedmann's equations obtained assum- 
ing a Robertson- Walker (RW) cosmological model and 
a perfect fluid moving along the geodesic congruence fol- 
lowed by the fundamental observers. It is well-known 
that stress-energy tensors corresponding to standard flu- 
ids cannot be responsible for the accelerated expansion. 
From a more general point of view, the Friedmann's equa- 
tion involving the acceleration of the scale factor results 
from the Raychaudhuri's equation assuming GR. This 
equation provides the expansion rate of a congruence of 
timelike or null geodesies (see [T5HPT] and recent review 
|18j). The Raychaudhuri equation plays an important 
role in the demonstration of the singularities theorems 
proved by Hawking and Penrose [13] . It is usually inter- 
preted that the contribution of space-time geometry to 
this equation represents the attractive (or non-attractive) 
character of gravity. An analysis of this contribution [TH] 
showed its geometrical interpretation as the mean curva- 
ture [2U] in the direction of the congruence. Besides, it 
can easily be verified that for a RW cosmological model 
with a negative deceleration parameter this contribution 
is positive, i.e. the mean curvature in the direction of 
the fundamental congruence turns out to be positive [19] . 
Hence, the attractive character of gravity vanishes. From 
this analysis, it is clear that the accelerated expansion of 
the Universe may be in conflict with the attractive char- 
acter of gravity. 



In GR, the attractive character of gravity is assured 
by assuming the usual energy conditions |13[ 1 14] . There- 
fore, a positive contribution to the Raychaudhuri equa- 
tion from space-time geometry is not attainable in GR 
provided that these energy conditions hold. Neverthe- 
less, this does not need to be the case in the context 
of modified gravity theories. In these theories, even if 
the usual energy conditions are assumed, a positive con- 
tribution to the Raychaudhuri equation from space-time 
geometry may be obtained. Moreover, an upper bound 
to the contribution of space-time geometry can be pro- 
vided both in terms of the gravitational model and the 
metric under consideration. Using this upper bound and 
assuming the usual energy conditions, throughout this 
investigation we shall derive restrictions to f(R) models 
in order to constrain their cosmological viability. 

Energy conditions have been widely studied in the lit- 
erature for different modified gravity theories. The au- 
thors of [2T] generalized energy conditions for a perfect 
fluid in f(R) theories by analogy with GR. In [52], the 
extended energy conditions of f(R) are used to derive 
energy conditions in Brans-Dicke theories with a van- 
ishing kinetic term in the Lagrangian using the equiva- 
lence between both theories. The energy conditions have 
also been studied for f(R) theories with a non-minimal 
coupling to matter [53]. In [53], the same procedure 
is applied to Gauss-Bonnet theories. In [55] and [55] . 
the authors considered Gauss-Bonnet theories with non- 
minimal coupling to matter and derive the correspond- 
ing energy conditions. All the aforementioned references 
followed the formalism first developed in [5TJ. This gen- 
eralization of the energy conditions, as the authors of 
[2T] themselves first acknowledged remains doubtful since 
there is no natural motivation but only an analogy with 
GR. This extension is in fact only motivated when the 
new terms appearing in the field equations are identified 
with physical fields. Nevertheless, these new terms may 
be understood as possessing only a geometrical meaning. 
Thus, there is no reason to assume any energy conditions 
on these terms. 

Moreover, if these new energy conditions are assumed, 
the mean curvature in every timclikc direction is non- 
positive by construction and, as already mentioned, for 
a RW space-time experiencing an accelerated expansion 
implies necessarily a positive mean curvature in the di- 
rection of the fundamental congruence. In particular, the 
stress-energy tensor associated with a cosmological con- 
stant A does not satisfy the usual energy conditions. This 
fact shows by itself the limitations of previous investiga- 
tions assumptions in the most trivial modified gravity 
Lagrangian beyond pure GR. 

This paper is organized as follows: First, in Section \H\ 
we summarize the standard energy conditions as consid- 
ered in GR and then the role of the Raychaudhuri equa- 
tion in the theorems of singularities is briefly discussed. 
Also, the procedure to be developed in the following sec- 
tion is also sketched. Section UTTl is devoted to introduce 
the field equations for f(R) theories in the metric for- 



malism as well as the commonly assumed conditions for 
f(R) models to be cosmologically viable. Then, in Sec- 
tion |IV] we assume the energy conditions in the frame- 
work of f(R) theories and present the inequalities that 
are obtained. We shall proceed by studying configura- 
tions of constant scalar curvature in Section [V] These 
configurations will enable us to impose some constraints 
on the parameters of several relevant f(R) models in or- 
der to get a positive contribution to the Raychaudhuri 
equation. Finally, we conclude our analysis by present- 
ing our conclusions. 

Throughout this study, we use a metric signature 
(—,+,+,+) and our definition of the Riemann tensor 
is: 

Rabc 1 = dbT ac — d a r^ c + r^ c T eb — r^,rf a , (i) 

R ac = R abc b holds for the Ricci tensor and R = R a a is the 
Ricci scalar. With this convention the usual Einstein's 
equations yield: 

Rab — n^-9ab = J— Tab ■ (2) 

From now on, we shall adopt c = G = 1. Furthermore, 
the stress-energy tensor corresponding to a perfect fluid 
is: 

Tab = p 6,6 + p (g a b + 6,6) . (3) 

II. ENERGY CONDITIONS IN GENERAL 
RELATIVITY 

The Raychaudhuri equation for timelike geodesies can 
be expressed as [TU [J5] 

~ = - \0 2 - Oab^ + U ab U: ab - RabC^ , (4) 

dr 3 

where 9, a a b and to a b are respectively the expansion, shear 
and twist of the congruence of timelike geodesies gener- 
ated by the tangent vector field £° and r is an affine 
parameter. On the other hand, the analogous equation 
for null geodesies becomes [HI ITS] 

(\B 1 

^ = -2 2 - a ab v ab + Lo ab u ab ~ R ab k a k b , (5) 

where k a is the tangent vector field to a congruence of 
null geodesies and A is an affine parameter. Let us recall 
that Q and ^ are geometrical identities, thus they hold 
independently of the gravitational theory assumed. 

In this investigation we are interested in the contribu- 
tion of space-time geometry to the previous equations, 
i.e., —R a b£, a £, b and —R a bk a k b . If a particular form for 
the metric tensor is assumed a priori, the Ricci tensor 
can be determined and then those contributions can be 
directly studied as it is the case for a RW metric [19] 
without considering any underlying gravitational theory. 
However, the metric tensor is in general unknown from 
the beginning, thus an expression for the Ricci tensor is 
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not at our disposal. In the latter scenario, the problem 
can be nonetheless tackled by using the field equations to 
obtain information about —R ab S, a £, b an d ~R a bk a k b . This 
procedure leads to considering conditions to be imposed 
on the stress-energy tensor T ab , the so-called energy con- 
ditions. Let us thus revise these conditions and some 
of their implications in GR both for timelike and null 
vectors. 



Timelike vectors 

The energy density of matter as measured by an ob- 
server with velocity £ a , is T ab £ a £ b . It is reasonable that 
this density would be non-negative. This requirement is 
known as the weak energy condition (WEC) 



T ab ^e > 



WEC. 



(6) 



Moreover, the dominant energy condition (DEC) ensures 
that the speed of the flux of energy is less than the speed 
of light, yielding 



T ab CT bc ( c < 



DEC. 



(7) 



which expresses that the flux of energy, i.e. — T ab £ a , is 
a timelike vector where the minus sign appears because 
we have chosen signature (—,+,+,+). Furthermore, we 
are mainly interested in the expression R ab £, a £, b - Using 
the usual Einstein's equations ^ and the subsequent 
relation between the Ricci scalar and the trace of the 
stress-energy tensor, i.e. R — — 8ttT, we obtain 



Rab^C = 8TT [ T a b ~ ^T 9ab ) 



a /-b 



8tt I T ab et b + =:T 



1. 



(8) 



It is customary to assume the positive sign of the r.h.s. 
of this equation since a distribution of standard matter 
would not result in a stress-energy tensor with pressure 
so large and negative as to make this member negative. 
This statement can be understood after replacing expres- 
sion ([3| in Hence, stress-energy tensors for standard 
matter fluids satisfy the so-called strong energy condition 
(SEC) 



TabCe > -\t, 



SEC. 



(9) 



It is known that both dust matter and radiation satisfy 
the SEC. For a discussion about cases where this condi- 
tion does not hold see |13j . In particular, a stress-energy 
tensor corresponding to a cosmological constant A fluid 
does not fulfill the SEC. We will discuss this case at the 
end of the section to avoid losing continuity in the dis- 
cussion. Therefore, the SEC requires 



iUCf >o, 



(10) 



which may be interpreted, because of asserting a non- 
positive contribution to Raychaudhuri equation, as a 



manifestation of the attractive character of gravity. It 
follows that the mean curvature [HI [20] in every timelike 
direction defined by 



-Rabi a i b 



(11) 



is negative or zero in GR provided that the SEC is as- 
sumed. 

The usefulness of the Raychaudhuri equation in the 
singularity theorems is based upon the following result: 
if one chooses a congruence of timclike geodesies whose 
tangent vector field is locally hypersurface-orthogonal, 
then uj ab — for all the congruence (as a consequence of 
Frobenius' theorem [Tl]) is obtained. The term a ab a ab 
is non-negative and whenever R ab £ a t; b > is assumed, 
then 



d(9 1 

d7 + 3 



9 2 < 0. 



which implies 



- 1 (r)>0o 1 + ^- 



(12) 



(13) 



This inequality tells us that a congruence initially con- 
verging (9q < 0) will converge until zero size in a finite 
time r < 3/|6> |, or in a reversed sense, if the congruence 
is initially diverging 9q > it was focused until zero size 
in the past. This result is important in proving the the- 
orems concerning singularities of Hawking and Penrose 
[131 H"4"] . Very often it is claimed that these theorems 
require energy conditions to hold, since for instance the 
SEC as we have just seen implies R a b£, a £, b > 0. However, 
these theorems are essentially mathematical theorems in- 
dependent of the gravitational theory. Energy conditions 
are necessary for these theorems to hold if and only if 
GR is assumed. Otherwise, the requirement would be 
Rab£, a £, b > for every non-spacelike vector. In Section 
IV[ we shall assume the usual energy conditions in the 
framework of f(R) theories but the sign of Rab£, a £, b will 
remain in principle undetermined. 



Null vectors 

Let us now consider a congruence of null geodesies. 
Just by replacing £ a — > k a in equation pi), one gets 

R ab k a k b = 8tt (r ab - \Tg a ^J k a k b = 8nT ab k a k b .(14) 

Hence, the so-called null energy condition (NEC) 

T ab k a k b > NEC , (15) 

implies that R ab k a k b will be non-negative. NEC is ful- 
filled by continuity if the SEC is assumed, but it is also 
fulfilled by imposing the WEC. Then, the assumptions 
for a congruence of null geodesies to focus are weaker 
than those for a congruence of timelike geodesies sketched 
above. Reasoning in the same way as before, if a congru- 
ence of null geodesies is initially converging 9q < it will 
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converge until zero size in a finite time r < 2/|0q| or in 
the reversed sense. 

The convergence of timelike and null geodesies in a 
finite time - in the future or in the past - is ensured 
in GR under the assumptions of the SEC and the NEC 
respectively. This result is usually known as the geodesic 
focusing theorem [131 1 14] . 



Beyond General Relativity 

If one considers a congruence whose tangent vector 
field is locally hypersurface-orthogonal, this means that, 
on the one hand LJ a t, — and, on the other hand, the r.h.s. 
of the Raychaudhuri equation for timelike geodesies Q, 
has only non-positive contributions of the parameters of 
the congruence (— 9 2 /3 and —a a ba ab ). If the observed 
acceleration of the Universe needs to be explained, a pos- 
itive contribution of the space-time given by Ai^a > is 
required at least for some directions £ a . This is unfeasible 
in GR if the SEC is satisfied as can be seen from (10 1. 



This opens two ways of circumventing the unavoidable 
attractive character of gravity in GR: either to suppose 
that the SEC is not satisfied or to modify the Einstein's 
equations. An analogous discussion may be done for null 
geodesies replacing the SEC by the NEC. 

In this work we shall consider the second scenario: the 
SEC will be located in a privileged place with respect 
to the Einstein's equations. We are thus assuming that 
standard matter satisfies the SEC and that the possibil- 
ity of a positive contribution to Raychaudhuri equation 
through the mean curvature M^a must be obtained from 
the f(R) modified field equations. By proceeding in this 
way, an inequality involving R a b^ a £, b and terms depend- 
ing on the gravitational theory under consideration will 
provide us an upper bound to the contribution of space- 
time geometry Ai^ . This bound will allow us to derive 
some restrictions on the f(R) models in order to get Ai^ 
positive, or equivalently R a b£, a £, b negative. 

As it was mentioned above, the cosmological constant 
stress-energy tensor does not satisfy the SEC and should 
not be regarded as an stress-energy term but as a partic- 
ular f(R) = -2 A model. 



III. f(R) THEORIES 
Let us consider the total action 



5' 



Sr, 



(16) 



i.e. a gravitational action plus a matter action term that 
includes all matter fields. The modification of f(R) the- 
ories to GR consists in assuming S grav of the form 



Sqrav , 

167T 



l - J (R + f(R)) V[i]d 4 x. 



where f(R) is an arbitrary function of R and g is the 
determinant of the metric. The variation of ( 16 1 with re- 



spect to the metric tensor and using (17 1 and (19 1 yields 



Rab{l + f{R))-\gab (R + f(R)) 



(V a V b - g ab n) f(R) = 8^ T ab . (18) 



where stress-energy tensor is defined as 

2 & Scatter 



Sg a b 



(19) 



Equations ( 18 ) are obtained in the metric formalism, i.e. 



the connection is assumed to be Levi-Civita connection. 
f(R) theories have been proved to be able to reproduce 
the cosmological history from inflation to the actual ac- 
celerated expansion era. For instance, it has been showed 
that the current evolution of the Universe can be repro- 
duced with certain f(R) functions [27] . 



(17) 



Let us remark here that equations ( 18 ) are fourth or- 
der differential equations. This is the reason of some 
strong instabilities that arise for certain f(R) models, 
such as the Dolgov-Kawasaki instability for the model 
f(R) = — [i A /R [28]. Moreover, there is a general insta- 
bility known as Ostrograski instability associated with 
Lagrangian that contains non-linear second derivatives 
terms. However, it has been proved that one can avoid 
Ostrograski instabilities in f(R) theories (cf. [29 ). The 
Cauchy problem in f(R) theories has also been consid- 
ered [30] where the authors concluded that the problem 
is well-posed in the metric formalism. Further details as 
well as local and cosmological tests for f(R) theories can 
be seen in [TTJ [T3] among others. 



Viability conditions of f(R) theories 

Some constraints are usually imposed on the f(R) 
functions in order to provide consistent theories of grav- 
ity. We would like to highlight two clear results that we 
will take into account in our study: 

1.1 + f'(R) > 0. This condition is imposed in order 
to ensure a positive effective gravitational constant 
G e ff = G/(l + f'(R)). It means that the main 
part of the contribution to the Einstein's equations 
conserves the sign [31]. This condition also guar- 
antees the non-tachyonic character of the standard 
graviton. 

2. f"(R) > 0. It ensures a stable gravitational stage. 
It is directly related to the presence of a positive 
mass in a high curvature regime for the scalar mode 
associated with this type of theories [32] . 

We will assumed both conditions for our discussion. 
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IV. ENERGY CONDITIONS IN f(R) THEORIES 

In this section we are interested in analogous equations 
to ([s]) and (14 1 when extended to f(R) theories. Thus, 
by imposing the usual energy conditions, inequalities in- 
volving the terms R a b£ a £ b (or R a bk a k b ) and f(R) extra 
geometrical terms will be obtained. These inequalities 
will provide us an upper bound for the contribution of 
space-time geometry to Raychaudhuri equation for time- 
like geodesies Q and for null geodesies (JjjJ. 

The usual approach in literature consisted of defining 
an effective stress-energy tensor by analogy with GR that 
includes the new geometrical terms in order to obtain an 
expression for R a b£ a £ b - Then, analogous energy condi- 
tions on this effective tensor were imposed. By this pro- 
cedure, a negative or zero contribution to Raychaudhuri 
equation from the term M.^ = —R a b£, a £, b for every time- 
like direction £° is obtained whenever these analogous 
energy conditions hold. Nevertheless, as we have already 
mentioned, M^a > is satisfied for almost all timelike 
directions in a RW cosmological model with the present 
value of the deceleration parameter q [TJJ]. Therefore, 
if those extended energy conditions hold, the present ac- 
celerated expansion of the Universe cannot be accommo- 
dated. 



Inequalities derivation 



Let us first take the trace of equation ( 18 ) that can be 
recast as 



R = 



-8nT-2f(R)+3Df'(R) 
1 - f'(R) 



(20) 



Therefore, from (181 together with (20) one gets 



Rab (1 + f'(R)) - \g ab (Rf'(R) - f(R) + 3D/'(JZ)) 



- (V a V b - g ab U) f(R) = 8ir \T ab - -Tg ab j . (21) 

Contracting the last equation with £ a £ b , where £° is a 
normalized timelike vector, £°£ a = — 1, we get 

RabCe (i + f'(R)) - (V£ h v a v b - ±n\ f'(R) 

+ l - (Rf'(R) - f(R)) = Stt (r ah Ci h + \t} • (22) 

If we consider a null vector k a instead of a timelike vector, 
then multiplying |2ll) by k a k b the result becomes 



R ab k a k b (1 + f(R)) - k a k b V a V b f(R) = 8TrT ab k a k b . 

(23) 

At this stage, let us impose the SEC and the NEC to 
the standard cosmological fluids in the expressions (22) 



and ( 23 ) respectively. After some manipulations, they 
1 



become 

RabCi" > 



2(1 + /'(A)) 

x [f(R) - Rf'(R) + (2C£ b V a V b - □) f(R)] (24) 



and 



R nb k a k b > 



l + f'(R) kakbWaWbf,{R) ' (25) 



where in both expressions 1 + f'{R) was assumed to be 
positive in order to guarantee G e ff = G/(l + f'(R)) > 0. 
Let us remind that our conclusions will be based upon 
this requirement. As a result of (24) and (25), one can 
conclude that although the SEC (NEC) have been as- 
sumed, in f(R) theories the sign of R a b£, a £ b (R a bk a k b ) 
cannot be determined a priori. Thus, if a certain model 



f(R) renders negative the right-hand side (r.h.s.) of (24) 
or (25), some freedom remains for RabCZ b or R ab k a k^to 



be negative which may be interpreted as a repulsive force. 
Reminding now the definition (11), the inequalities (24) 



and ( 25 ) can be cast in the following way 



M*a = -R ab CC < 



-1 



2(1 + /'(£)) 

[f(R) - Rf'(R) + (2re b V a V b - □) f'(R)} , (26) 



R ab k a k b < 



-1 



i + f'(R) 



k a k b V a V b f'(R) ; (27) 



that provide upper bounds to the contribution of space- 
time geometry to the Raychaudhuri equation for timelike 
and null geodesies respectively. 



Let us stress that in vacuum the inequalities (24) and 



( 25 1 - or equivalently ( 26 1 and ( 27 ) - get saturated since 



the SEC/NEC energy conditions are trivially saturated. 
Consequently, if a f(R) model renders the r.h.s. of (26) 
and ( 27 ) positive in vacuum scenario, then a positive con- 



tribution to the Raychaudhuri equation is automatically 
obtained for timelike geodesies and null geodesies respec- 
tively. 

Let us focus on a short example. The Einstein's equa- 
tions with a cosmological constant A become 



Rab 



1 



-Rg a b + Agaf, = 87TT ah 



(28) 



As we have previously noted, GR with a cosmological 
constant is equivalent to take f(R) = — 2A. In this trivial 
case, the inequalitites ( 26 ) and p7|) become 



Me* < A 



-R ab k a k b < 0. 



(29) 



The first inequality tells us that in the case of timelike 
geodesies a positive contribution to the Raychuadhuri 
equation from space-time geometry M^a > is possible 
provided that A > 0, which corresponds to the correct 
sign of A to provide cosmic acceleration. 

Hence, it may be thought that a simple criterion to 
decide when a f(R) model is able to render accelerated 
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expansion of the Universe has been obtained. Nonethe- satisfies 



less, expressions ( 24 1 and ( 25 ) need to be evaluated at the 



solutions for ( 18 ) so the problem remains cumbersome. 



However, such a problem is absent when space-time con- 
figurations with constant scalar curvature Rq are consid- 
ered. 



Constant scalar curvature solutions 

Space-times both in vacuum and GR cosmological sce- 
narios when studied at late times, with both radiation 
and dust being negligible with regard to a cosmological 
constant, are maximally symmetric, i.e. they possess a 
constant Gaussian curvature Kq. This implies a con- 
stant scalar curvature Rq — 12Kq (but the reverse is not 
generally true) . For this reason, it may be expected that 
solutions of constant scalar curvatures will be recovered 
at late times by physically viable f(R) models. 

Since the covariant derivatives of f'(R) in solutions of 
constant scalar curvature (R = Rq) are zero, expressions 
([24l and (l25|) result respectively in 



Rab^e > 



/(go) - Rof'(R ) 
2(l + /'(i? )) 



and 



R nh k a k b > 0. 



(30) 



(31) 



A remarkable result follows from the inequality (31 1: In 



f{R) theories, the condition for the null geodesic focusing 
theorem to hold, namely R ao k a k b > 0, is satisfied in 
space-times of constant scalar curvature provided that 
the NEC is assumed as given by ( 15 ). It is worth noticing 



that this result does not depend upon the sign of Rq nor 
upon the f(R) model under consideration. Moreover, 
the NEC is only assumed in the standard stress-energy 
tensor for matter, not in the effective one usually defined 
after gathering all the new terms of the modified Einstein 
equations. Since the holographic principle 33 makes use 
of the null geodesic focusing theorem in order to ensure 
that light-sheets will eventually end, the previous result is 
of extraordinary importance when studying this principle 
in f(R) theories. 

In the rest of this investigation we shall focus on time- 
like geodesies. Therefore, the r.h.s. of (30) must be 
negative in order to allow R a b£, a £, b < or cquivalently 
A4^a > 0. Thus, Mta be bounded from above. Hence 
we impose 



f(Ro) - Rof'(Ro) 



<0, 



2(l + /'(i? )) 
and provided that G e // > 0, we get 

f(Ro)-Rof'(R o )<0. 



(32) 



(33) 



If we now consider the equation |20[ ) in vacuum (T = 0) 
for constant scalar curvature solutions, the value of Rq 



Ro 



-2/(flp) 

i - f'(Ro) 



(34) 



which is an algebraic equation relating Rq with the pa- 
rameters of the f(R) model under study. Although in 
general this equation cannot be solved analytically, there 
exist some f(R) models for which a closed solution de- 
pending upon the parameters of the model can be found. 
Using the equation |34|) in ( 33 1 one gets 



/(-Ro 



1 - f'(Ro) 



<o, 



that together with ( 34 ) , implies 



Ro > 0. 



(35) 



(36) 



Hence, a positive contribution to the Raychauhduri equa- 
tion from the space-time geometry M^a for every time- 
like direction is obtained provided that R > 0. This 
condition will constrain the parameters of different f(R) 
models as will be seen in the next section. 

In fact, there exists another straightforward way of get- 
ting ( 35 1 as follows: The solutions of ( 18 ) in vacuum with 



constant scalar curvature imply 



1 



R ab (1 + f'(Ro)) ~ ^ 9ab (Ro + f(Ro)) = (37) 



and consequently 



Rab 



1Rq_ 
2 14 



f(Ro) 



f'(Ro) 



9ab 



Rq 



~9ab ■ 



(38) 



where (34) has been used in the last equality. It means 



that the allowed f(R) space-times with constant scalar 
curvature in vacuum are Einstein spaces [20 . Thus, if 
a negative value of R a b£, a £, b is required in order to have 
M^a > the condition R > needs to be accomplished. 
However, one must keep in mind that the inequality (33) 



is more general since it allows us to have > even 
for non- vacuum scenarios. 

If one consider a maximally symmetric space-time, the 
condition ( 36 1 implies that in order to get a positive 



contribution to the Raychadhuri equation the space-time 
must be de Sitter (Rq > 0). The constraints on the pa- 
rameters of different f(R) models explored in the follow- 
ing guarantee that such a space-time exists as a solution 
of the modified Einstein's equations (18). 



V. f(R) MODELS 

Let us now study some f(R) models in vacuum to il- 
lustrate the previous results: 

• Model I f(R) = a\R\ P 

This model encompasses a wide variety of propos- 
als available in the literature. For instance, the 



6 





IIP 











a 



Figure 1. (a, /?) plane for Model I in vacuum: The pairs (a, (5) 
which lying in the blue region are those that satisfy a(/3 — 2) > 
and consequently Rq > 0. The meshed zone fulfills f"(Ro) > 
which is a stability condition commented in Section |III| The red 
zone does not satisfy the condition G e jf = G/(l + /'(Ho)) > 
and therefore the inequality | |35| l is not valid there as we have 
already commented in the previous section. The parameters that 
provide a positive contribution from the space-time geometry, i.e. 
JVi^a = —R a b^ a i b > 0, to Raychaudhuri equation for congruences 
of timelike geodesies are those of the blue zone excluding the red 
zone about which no statement can be done with our discussion. 



case with (3 = 2 has been considered by Starobin- 
sky in relation with the inflationary epoch in [34] 
and more recently in [35] as a candidate to take 
account of the dark matter. The case (3 < 1 and 
in particular (3 = — 1 was proposed in |3BJ as a 
possible mechanism to provide cosmological accel- 
eration, although it is currently ruled out. For this 
model, the curvature scalar can be obtained from 
(34 1 yielding 



Ro 



± 



±1 



a(fi-2) 



(39) 



where the sign depends upon the sign of Rq because 
of the derivative of the absolute value (plus signs for 
Rq > and minus signs for Rq < 0) . There is also a 
trivial solution with Rq = which is of no interest 
for our discussion. Since we are interested in Rq > 
we take the expression with the plus signs. In 
order to a positive constant scalar curvature exist 
the parameters must obey 



a{P-2) > 0. 



(40) 



Note that for (3 = 2 the only constant scalar curva- 
ture in vacuum is Rq = 0. 



Regions where the conditions a((3 — 2) > and 
/"(-Ro) > hold are plotted in Figure[l] The region 
where G eff = G/(l + f(Ro)) > does not hold is 
also represented. Since this viability condition of 
f(R) has been assumed in deriving the inequality 
( p5| , our discussion is not valid for the values of the 
parameters falling in that region. Let us stress that 
the conditions /"(Rq) > and G e // > are eval- 
uated in the corresponding value of the constant 
scalar curvature Rq which depends on the parame- 
ters. It means that these conditions will be satisfied 
by the parameters that fall in the corresponding re- 
gions in the case of constant scalar curvature solu- 
tions, not for every solution of the equation (18 1. 



This consideration remains valid for all the models 
studied in this section. In Figure [l] one can also see 
that for (3 > there are regions where both require- 
ments, namely Rq > and /"(Rq) > 0, hold. For 
a < 0, (3 must be restricted to the interval (0,1); 
on the contrary, for a > 0, it must be (3 > 2; in 
order both requirements to hold. 

Model II f(R) = R a exp((3/R) - R 

This model was discussed for a = 1 in [37] and 
more recently in [35]. By proceeding as for the 
previous model, one gets from equation (34) 



Rq — 



(41) 



Thus, in order to ensure a positive contribution to 
the Raychaudhuri equation from the space-time ge- 
ometry, i.e. = —R a b£, a £, b > one must impose 







> 0. 



(42) 



Note that for a = 2, the only constant scalar cur- 
vature vacuum solution is Rq = 0. The same condi- 
tions as for the previous model are plotted in Figure 
[2] Let us remind that the conditions G e f / > and 
j"(Ro) = are evaluated in Rq which depends on 
the parameters of the model. 

Let us stress that for this case the value of f"(Ro) 
depends on the sign of Rq. Since we are interested 
in the region where Rq > holds, the condition 
/"(-Ro) — is plotted after assuming R > 0. For 
this model, in the zone where G e // > holds, all 
the parameters that provide a positive scalar cur- 
vature Rq > also fulfill /"(-Ro) > 0. 



Model III f(R) = R [log(aR)f - R 

This model was also considered in [371 1 
case, equation renders (34) 



Rq 



■ exp(/3) . 



In this 



(43) 
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Figure 2. (a, /3) plane for Model II in vacuum: For the parameters 
in the blue zone a positive constant scalar curvature exist. The 
condition /"(Ro) > is fulfilled in the meshed zone. Let us remark 
that in this model the condition f"(Ro) > depends on the sign 
of Rq and the value of a. For simplicity, the condition /"(Ro) > 
is plotted after assuming Ro > since it is the case in which we 
are interested. Moreover, if Ro > the condition /"(-Ro) > 
does not depend on the value of a. In the red region the condition 
G e ff = G/(l + /'(-Ro)) > does not hold and consequently the 
inequality ||35| does not apply there. 



Therefore, for this model the condition guarantee- 
ing both Rq > and a positive contribution to 
Raychaudhuri equation for timelike geodesies from 
the space-time geometry is a > 0. 

The same conditions considered for the previous 
models are plotted in Figure [3] For this case, there 
exists also a region where both conditions Rq > 
and /"(Rq) > are satisfied, namely for a > and 
> 1/2. 



Model IV f(R) 



This model has been proposed in as cosmo- 
logical viable attracting much attention in the last 
years. In order to illustrate our procedure, let us 
consider the particular case with n = 1 



f(R) = 



aR 

i + Jr 



(44) 



where a trivial redefinition of the parameters has 
been performed. The constant scalar curvature of 
this model in vacuum becomes 



R± 



P 



1 ± y/a(a-l) 



(45) 



Figure 3. (o, 0) plane for Model III in vacuum: As for the previous 
figures, Ro > is obtained with parameters in the blue zone. In 
the meshed zone /"(R ) > 0. G eff = G/(l + /'(Ro)) > is not 
fulfilled in the red zone and thus the inequality | |35| l is not valid 
there. Let us remark that the condition /"(Ro) > depends on 
the sign of 0. As we are interested in region where G e ff > holds, 
it is plot the region where /"(Rq) > assuming > 0. 



It follows that only for a(a — l) > constant scalar 
curvature solutions exist. Therefore, imposing the 
constraint R± > we get 



1 y/a(a- 1) 



P 



> 0. 



(46) 



Analogous plots to the previous models are shown 
in Figure [4j For this model, two different figures 
are shown since there are two possible values of Rq, 
namely R±. For the case the regions J?_ > 
and f'(R-) > do not overlap. On the other hand, 
for the case all the parameters pairs providing 
a positive constant scalar curvature R + > also 
fulfill f"(R + ) > 0. Furthermore, for R + case, the 
condition G e // is always satisfied. 



VI. CONCLUSIONS 

In this investigation we have studied the contribution 
of space-time geometry to the Raychaudhuri equation for 
timclike and null geodesies. In General Relativity with- 
out a cosmological constant, once the usual energy con- 
ditions are assumed it is not possible to obtain a pos- 
itive contribution from the space-time geometry to the 
Raychaudhuri equation for timelike and null geodesies. 
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Figure 4. (a, /J) plane for Model IV in vacuum: i?+ is considered 
at the upper panel whereas the lower panel considers R-. The 
blue zone represents the region where R+ > (R— > 0). In the 
meshed zone the condition f"(R+) > (f'(R-) > 0) holds. The 
red zone parameters can not be considered in our discussion since 
G ef f = G/(l + f'(R+)) > (G eff = G/(l + f'(R-)) > 0) docs 
not hold there. Finally, the parameters that fall in the grey zone 
do not fulfill a(a — 1) > which is a necessary condition in this 
model in order to have a solution with constant scalar curvature. 



Relativity are assumed, the fact of getting a positive con- 
tribution to the Raychaudhuri equation from space-time 
is allowed. 

We have derived two inequalities that bound from 
above this contribution for congruences of both timelike 
and null geodesies. In order to allow a positive contribu- 
tion to the Raychaudhuri equation, these upper bounds 
must be positive. The limitation with the obtained in- 
equalities is that in general in order to extract some infor- 
mation, a metric solution of the modified Einstein equa- 
tions must be used. Nevertheless, in the cosmological 
relevant case of constant scalar curvature i?o solutions, 
such a knowledge is not required. Under this assump- 
tion, it was obtained that R a i,k a k b > where k a is a null 
vector. This is the condition needed for the null geodesic 
focusing theorem to hold. Thus, this theorem remains 
valid in the context of f(R) theories for space-times with 
constant scalar curvature, regardless the f(R) model con- 
sidered and the matter content (provided that the Null 
Energy Condition is fulfilled). This result acquires a re- 
markable importance when dealing with the holographic 
principle in the context of f(R) theories. Since if the null 
focusing theorem holds for this scenario in fourth order 
gravity theories, the light-sheets will eventually end. 

Finally, in vacuum scenarios for space-times of con- 
stant scalar curvature Rq, we derived constraints on the 
parameters of paradigmatic f(R) models guaranteeing a 
positive contribution to the Raychaudhuri equation. We 
conclude that for the models under consideration, there 
exist parameters values that allow the desired contribu- 
tion while satisfying the condition /"(Rq) > 0, which 
ensures the stability of the solutions. 

In this investigation, usual energy conditions were 
solely imposed upon the cosmological standard fluids 
whereas the extra f(R) terms were considered as geomet- 
rical terms. This approach renders a more realistic anal- 
ysis of the viability of f(R) models and has been studied 
here for the first time. Furthermore, the sketched proce- 
dure developed constitutes a straightforward and system- 
atic approach to decide whether a particular f(R) model 
could generate cosmological acceleration. For this pur- 
pose, we have paid particular attention to constant cur- 
vature solutions. With the tools presented in this inves- 
tigation, analysis can be extended to more complicated 
cosmological scenarios and other alternative gravity the- 
ories beyond the Concordance model |40j . 
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